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SUMMARY 


Background stratospheric aerosols typically yield scattering ratios R of 
the order of 1.10 at 20 km for a ruby lidar system. Background aerosol amounts 
inferred from lidar backscatter measurements are affected by both random and 
systematic uncertainties. In this report, the calibration and data-analysis 
procedures for the NASA Langley 48-inch lidar system are discussed. These 
procedures have been supported by obtaining concurrent dustsonde and lidar 
profiles of stratospheric aerosols during a comparative experiment conducted in 
San Angelo, Texas, in May 1974. A comparison of profiles obtained by these two 
techniques is within the accuracy of the two measurements below 23 km. At 
higher altitudes the comparison is improved when an altitude-dependent aerosol 
size distribution is considered. Lidar data from the Texas experiment yielded 
an average random uncertainty in R of 2.5 percent, with uncertainties in 
molecular density being most important. The systematic uncertainty in R due 
to the lidar profile normalization procedure is estimated from aerosol model 
calculations to have been 5 percent. For the Texas experiment the aerosol 
component of the atmospheric backscatter at 20 km contains random and system- 
atic uncertainties of 35 and 63 percent, respectively. 


INTRODUCTION 

Various remote measurement techniques are being used to characterize the 
stratospheric-aerosol layer. One such technique is the use of lidar systems. 

A lidar measure of the aerosol concentration is the scattering ratio R, that 
is, the total atmospheric backscatter divided by the molecular backscatter. 

The quantity R - 1 is the fraction of the atmospheric backscatter ing that is 
due to aerosols. Efforts to monitor the seasonal and annual variability of 
stratospheric aerosols have been undertaken at a number of lidar installations. 
Fox et al. (ref. 1), Remsberg and Northam (ref. 2), and Russell et al. (ref. 3) 
have presented data on R - 1 from 1964 to 1974 to show the decreasing trend 
in aerosol loading at the 20 km altitude since the eruption of Mt. Agung in 
1963. (See fig. 1.) In figure 1 the vertical bars represent the range of 
scattering-ratio data over the time periods within, the hor-izontal bars. A 
background level of 1.20 to 1.05 in R prevailed from 1971 to 1974. Up to 
1974 the aerosol loading had either leveled off or was still declining. The 
occurrence of the Fuego volcanic event in October 1974 has made it difficult to 
know whether or not a steady-state stratospheric-aerosol level had been reached 
during the summer of 1974. Thus, the term "background aerosol amount" is still 
a relative one. The decay of the Fuego aerosol layer has been discussed by 
McCormick et al. (ref. 4); during 1977 the layer approached the pre-Fuego 
levels once again. Studies of trends in the background aerosol amount require 
data sets taken periodically at a number of geographic locations. 

The accuracy of the lidar measurements in figure 1 and the ability to 
intercompare them depend heavily on the individual system calibrations and 
data-analysis procedures. Northam et al. (ref. 5) have reported on a Wyoming 



comparative experiment (WCE) at Laramie, Wyoming, in 1972 using the Langley 
48-inch lidar system and the University of Wyoming dustsonde measurements. 

That particular experiment indicated agreement between the relative profile 
shapes obtained by the two methods, but further refinements in the lidar cal- 
ibration and data-analysis procedures were necessary before systematic error 
estimates could be made. Further calibrations of the Langley 48-inch lidar 
system were performed, and a Texas comparative experiment (TCE) was conducted 
on May 7 and 8, 1974, at San Angelo, Texas. This report describes the refine- 
ments in the data-reduction procedures employed in that experiment, and the 
various error sources are described in some detail. The results of the TCE 
are then discussed, including stratospheric-aerosol profiles as measured by a 
balloon-borne in situ particle counter and by backscattered light from the 
ground-based lidar operating at a wavelength of 0.694 ym. Aerosol extinction 
estimates derived from the lidar profiles are also compared with extinction 
data inferred from searchlight measurements of the aerosol layer. The implica- 
tions for the monitoring of long-term aerosol trends by lidar are reviewed. 

The authors wish to express their gratitude to Professor J. Rosen of 
the University of Wyoming for conducting the dustsonde launch at San Angelo, 
Texas, and for supplying his reduced data. M. Patrick McCormick and William H. 
Fuller, Jr., of the Langley Research Center participated in the coordination of 
the TCE, and W. Hunt, C. Bartusiak, and F. Diehl of Wyle Laboratories, Hampton, 
Virginia, assisted with the lidar measurements at both San Angelo, Texas, and - 
Hampton, Virginia. Dr. S. K. Poultney of Perkin Elmer Corporation, Norwalk, 
Connecticut (formerly of the Old Dominion University Research Foundation) and 
W. Hunt performed the photomultiplier-tube calibrations by simulating the lidar 
returns with light-emitting diode pulses. 


SYMBOLS 
E 
f 
k 
N 
ft 

P(180) backscatter ing phase function' 

q transmissivity of model atmosphere 

R scattering ratio 

r particle radius, ym 

T turbidity 

T(180) backscatter turbidity 


laser output energy per pulse, J 
backscatter ing function, m -1 -sr -1 
constant in equation (1) 
number density, m - 3- 
• refractive index 
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V voltage of photomultiplier, V 

Z altitude, m 

a aerosol size parameter 

3 total scatter function, m _1 

Y relative backscatter efficiency 

f| geometric cross section, m* 

X wavelength, ym 

v size distribution parameter 

a backscatter ing cross section, m^-sr -1 

fl solid angle, sr 

Subscripts: 

a aerosol 

d dustsonde 

m molecular 

Abbreviations: 

TCE Texas comparative experiment 

WCE Wyoming comparative experiment 


DESCRIPTION AND CALIBRATION OF LIDAR SYSTEM 

The Langley 48-inch lidar system (refs. 5 and 6) consists of a pulsed 
laser transmitter and a telescope receiver with parallel alignment of the 
optical axes. (See fig. 2.) The laser pulse of nearly monochromatic light 
is scattered by aerosols and molecules as it propagates through the atmo- 
sphere. A small percentages of this light is backscattered into the telescope 
which directs the light through optical filters onto one or more of the four 
photomultiplier tubes. The beamsplitter directs fractions of the total inci- 
dent light energy to the first three tubes, while the fourth tube is a photon 
counting detector. The detector package is triggered at the time of the laser 
firing, and each photomultiplier tube measures the backscattered photons or 
energy at appropriate subsequent time intervals. In this manner a range 
resolved measure of the total backscatter, due to both aerosols and molecules, 
is obtained. 
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Figure 2.- Langley 48-inch lidar system. 




The transmitter is a temperature-stabilized ruby laser, operating at 2.5 
to 3.0 J at 4 pulses per minute. The transmitter is equipped with a dou- 
bling crystal, allowing for two wavelength outputs (0.347 and 0.694 ym) . A 
neodymium-glass laser (1.06 and 0.53 ym) is available as an auxiliary trans- 
mitter. This report discusses only returns obtained at 0.694 ym. The receiver 
is an f/10 Cassegrain 1.32-m (48-in.) diameter telescope. The total system 
also includes a pulse counter, analog- to-digital converter, tape recorder, and 
computer subsystems. A calibrated photodiode samples a portion of the laser 
output energy and also prpvides a zero time pulse that initiates the system 
timing circuits. The photomultiplier outputs are amplified, analog-to- 
digitally converted at a rate of 10 MHz with 8-bit 2048-word capability, and 
placed in buffer memory. The computer accesses the buffer memory at a 1 MHz 
rate and stores the data on magnetic tape. 

In the WCE, complete profiles were obtained by combining three data segr 
ments which were taken at successive times. That technique imparted a potential 
error to the profile due to the temporal changes occurring in the atmospheric- 
aerosol returns from each altitude segment. To overcome this problem, a sequen- 
tially switched three-tube high-dynamic-range analog detector package was used 
each laser firing in the TCE to survey the 10 to 30 km altitube profile. Poten- 
tial nonlinearities in the detector tubes and amplifiers were detected by simu- 
lating the backscattered laser pulse with a light-emitting diode and by compar- 
ing the three channels. Linear tube response is a necessity for combining the 
signals measured with the three tubes and also for the subsequent normalization 
of the total profile to the expected molecular return at a given altitude. As 
pointed out by Pettifer (ref. 7) , signal- induced noise results in a residual 
nonlinearity that may have a significant effect on the final lidar profile. 

This effect is fully discussed for the Langley 48-inch lidar system by Hunt and 
Poultney in reference 8. Photomultiplier-tube response was checked for nonlin- 
earities both before and after , the TCE. A calibration of the 48-inch lidar 
data-acquisition system was also performed. Due to a faulty electrical align- 
ment, the analog- to-digital converter (ADC) response was nonlinear, but it was 
accounted for by calibration of the ADC system. 

Up to 10 background measurements were performed before and after each set 
of 50 to 100 laser firings in the TCE. Part of the background signal resulted 
from the ruby fluorescence return due to aerosols at low altitudes. This back- 
ground component occurs because the laser and telescope for the 48-inch system 
are separated by 1 m; this proximity resulted in overlapping of the transmitted 
and received images at Z = 1/2 km for a 4 mrad field of view. This low over- 
lap altitude was just where aerosols tend to have their highest concentration in 
the mixed layer and, thus, where background fluorescence was a problem. That 
part of the return was simulated by firing the laser at reduced Pockel cell 
voltages for which no lasing pulses were emitted. Those background returns 
varied slightly from one set of laser firings to another; and, since they were 
superposed on the stratospheric returns at altitudes where profile normalization 
was often attempted (Z = 28 to 31 km), adjacent sets of background returns were 
averaged and subtracted from each corresponding set of atmospheric returns. 
(Since the TCE, a rotating shutter has been installed to eliminate .the fluores- 
cence return.) Another source of background signal was the slowly varying "dc 
drift" due to the dc coupled amplifier and to low light levels being passed by 
the spectral filters. An uncertainty of ±1/2 count was also present for the 
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system itself because of an inability to record fractional counts. At the 
upper end of each of the three photomultiplier-tube profiles, the ratio of 
actual received signal to that measured as dc background became small. Any 
drift in the dc level was therefore important and was considered by occasion- 
ally comparing photomultiplier-tube output curves for any changes in gain. 


DATA ERROR ANALYSIS 

A typical lidar-return curve from the atmosphere for the TCE is given in 
figure 3. This curve represents range-corrected relative signal Z^v/E as a 
function of range Z, where V is the voltage of the return signal and E is 
the laser output energy. The molecular-return curve in figure 3 represents the 
calculated return signal which would be expected from just the molecules in the 
atmosphere at the time of the TCE. The difference between the two curves, 
then, is a measure of the aerosols or clouds that were present. Since the 



Figure 3.- Range-corrected lidar returns from TCE. 0330 to 0400 MDT; 100 shots. 
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atmospheric- signal profile V(Z) varies by over three orders of magnitude from 
4 to 30 km, three photomultiplier tubes were required to provide the dynamic 
range necessary to measure the backscattered signal accurately. Techniques for 
generating complete profiles (fig. 3) from those measurements are discussed in 
the following section. 


Background and Combination Effects 

At the time of the WCE (ref. 5) , the data- acquisition program limited the 
data return to 700 words, corresponding to a range increment of 10.5 km. For 
returns up to 30 km, it was necessary to record sets of three 10.5 km altitude 
segments of data at successive times by changing the initial delay on the 
photomultiplier-tube detector package. In order to obtain a single continuous 
profile, the altitude regions of overlap between sets of adjacent segments were 
used to combine those averaged segments with each other, and the result was fit 
to a polynomial. In the overlap regions, any pronounced features either due to 
fluorescence and detector background or to aerosol layer structure could have 
produced significant errors in the final combined profile. Real structural 
changes in the atmospheric returns, though infrequent, have been observed dur- 
ing a period of several hours. 

The sequence of data handling for the TCE proceeded as follows. A multi- 
plexer was used to switch-in the three photomultiplier tubes, and 2048 words 
were recorded to yield a total 30 km profile for each laser firing. Each seg- 
ment was corrected for variations in background effects and laser power. To 
determine background effects, 10 measurements were taken and averaged. Since 
signals from each of the three . photomultiplier tubes decrease with range, the 
ratios of signal to background approach unity and, therefore, the data become 
increasingly more sensitive to background subtraction (fig. 4) . As a result, a 
linear approximation was made to the second half of these returns to eliminate 
the effects of system noise fluctuations in the low-level background signal. 
Then, this approximation to the background was subtracted from the averaged 
set of stratospheric returns, and the resultant profile segments were range 
corrected. Due to the different amplifier gains and light levels for the three 
photomultiplier tubes, the data segments still had to be combined into one com- 
plete profile. To minimize uncertainties in atmospheric returns at the gain- 
switched altitudes, two sets of returns (primed and unprimed segments in fig. 5) 
were taken at about 15-min time intervals with combination points at different 
altitudes. Continuous portions of one profile overlapped the gain-switched 
altitudes of the other profile, and any signal perturbations were more easily 
discerned by comparing both profiles. The problem of rapidly (within 15 min) 
changing signal perturbations between adjacent sets of data still remained, 
however, as in the case of a thin cirrus feature moving across the field of 
view, visual inspection of the data at the time of the measurement identified 
those few cases, and profile data were obtained then with just a single photo- 
multiplier tube over a more limited altitude range. Once the data from the 
three photomultiplier tubes were combined, the profile was further averaged, 
this time over altitude. Typically, 30 altitude points were averaged to give 
a height resolution of 0.45 km. No other smoothing was applied to the data to 
achieve the final profiles. 
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Figure 4.- Background and lidar returns for three altitude ranges 



Figure 5.- Procedure for combining return signals. 



Random and Systematic Errors 


Northam et al. (ref. 5) have discussed the lidar equation and its subse- 
quent use in the calculation of scattering ratio R and aerosol backscatter ing 
function f a . The ratio R of the signal given by the lidar equation to the 
calculated molecular or clear-air return is given by 


Z 2 V(Z) f a< z ) 

R = k = 1 + (1) 

f m (Z)q 2 (Z) f m( z ) 


where q 2 is a two-way atmospheric attenuation factor, Z is altitude, V is 
return signal, and f m and f a are the molecular and aerosol backscattering 
functions, respectively. The quantity fm (Z)q 2 (Z)/Z 2 is proportional to the 
expected molecular lidar return. Since ,f m = o m N m (Z) , where N m is molecular 
number density and a m is the molecular backscattering cross section, the aero- 
sol portion of the scattering ratio f a (Z)/f m (Z) can be calculated. Since it 
is difficult to determine remotely whether any altitude region is free of aero- 
sols, the constant k that relates the expected molecular return to the total 
lidar signal return was evaluated on the basis of assuming only molecular scat- 
tering for altitudes between 28 and 31 km where aerosols often are at a minimum 
(refs. 3 and 5). This assumption forces R to be unity at those altitudes. 

Figure 6 shows the scattering-ratio profile derived from data in figure 3 
for an average of 100 lidar returns. This averaging procedure gives better than 
2-percent standard deviation in the scattering ratio for altitudes up to 30 km 
over 25 min. The molecular-return profile was calculated by using temperature 
and pressure data obtained at 0200 MDT at- San Angelo, Texas. Error bars in fig- 
ure 6 represent 95-percent confidence intervals for the random changes in the 
scattering-ratio data. These intervals do not increase uniformly with altitude 
because the signal-to-noise ratio for each photomultiplier tube was about the 
same. A bulge in the scattering-ratio profile occurs between 1 8 and 24 km, but 
its peak value is only 1.07 which is representative of relatively low aerosol 
content. Equation "(1) can be rearranged and used to calculate the aerosol 
backscattering function 


f a (Z) = (R - l)a m N m (Z) 


( 2 ) 


which is plotted in figure 7. Since f a - may be represented as the product 
of an effective aerosol backscatter cross section o a and particle concentra- 
tion N a (ref. 5) , profiles developed herein may be used to infer the aerosol 
number density as a function of altitude provided a a is known. The conver- 
sion of the lidar returns to a scattering-ratio profile and then to an aerosol 
backscattering-f unction profile was subject to several possible errors: mea- 

sured temperature profile, attenuation model, and normalization procedure. 

Each of these is discussed in turn. • - - 






The effect of using different molecular-density profiles for calculating a 
scattering-ratio profile is shown in figure 8 for another lidar data set taken 
on May 5, 1974. Radiosonde data were obtained at both San Angelo, Texas, and 
Midland, Texas {about 1 50 km to the northwest) about 3 hr after the lidar mea- 
surements. Due to differences in temperature and pressure from the two loca- 
tions, R was found to vary from 0 to 3 percent in this profile comparison, 
with the average being 1 percent. This uncertainty may, in fact, be due to the 
relative radiosonde density error of 1 percent (ref. 3) . Thus, uncertainties 
in air density are comparable to the 2-percent random signal error. Since the 
maximum contribution to R from aerosols was about 7 percent in the strato- 
sphere at that time, variations in molecular density were significant. 

Intervening atmospheric species may absorb a portion of the laser beam as 
well as scatter the incident light in all directions. Attenuation profiles for 
a particular wavelength are proportional to the number density of neutral mole- 
cules, aerosols, and ozone; the latter is especially important for wavelengths 
shorter than 0.32 Urn. Molecular scattering and ozone absorption remain nearly 
constant from year to year, but the aerosol or Mie scattering depends on the 
current particle size distribution and composition. 

Values for the two-way attenuation q^(Z) applied in figure 6 were 
approximated from Elterman's tabulation at 0.70 ym (ref. 9). No attempt was 
made to correct the Elterman attenuation data for the aerosol content being 
measured because simulations indicate that for aerosol amounts in the strato- 
sphere (R = 1.10) just before May 1974 an error in R of only 1 percent 
occurred if the aerosol attenuation was neglected altogether (ref. 3). 
Elterman's aerosol attenuation coefficients (ref. 10) derived from data taken 
after the volcanic eruption of Mt. Agung in 1963 are not applicable to an anal- 
ysis of background aerosol levels. 

The random uncertainties are summarized as a signal error of 2 percent, a 
molecular density error of 1 percent (range of 0 to 3 percent) , and an attenu- 
ation error of 1 percent. The root-mean- square value of these uncertainties 
equals 2.5 percent. 

In determining the magnitude of the scattering ratio with equation (1), 
one must assume that f a has a known minimum value (normally taken to be 
zero) at some altitude. Russell et al. (ref. 11) found that the region just 
above the tropopause at Laramie, Wyoming, is often aerosol free, but this char- 
acterization cannot be generalized for other locations because the occurrence 
of aerosol-free regions near the tropopause is effectively influenced by the 
tropospheric weather patterns and the position of the jet stream. Theoretical 
one-dimensional studies of processes affecting the background stratospheric- 
aerosol layer indicate that the major constituent, sulfuric acid aerosols, 
should be depleted because of "washout" near the tropopause and because of 
evaporation above about 30 km (ref. 12). In the present analysis the region 
between 28 and 31 km was selected as the normalization region. That choice is 
examined further in the following section. 
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TEXAS COMPARATIVE EXPERIMENT 


Lidar-Dustsonde Comparison 

The time line for the lidar measurements and for the dustsonde launch for 
the TCE is shown in table I. Over 800 laser firings were made on a schedule 
to include the dustsonde launch and measurements. The dustsonde, launched at 
0200 MDT, ascended through the lower stratosphere at approximately 0330 to 
0400 MDT in coincidence with the averaged lidar data in figures 6 and 7. The 
dustsonde measured the particle number density N<j directly using two detector 
channels, one (channel A) to discriminate particles larger than 0.15 ym in 
radius and the other (channel B) for particles larger than 0.25 ym. For each 
channel, 500 particle counts were registered before a data point was recorded; 
the counting error in N^ according to Poisson statistics was 4 percent. 
Hofmann et al. (ref. 13) give an uncertainty in N^ for channel A of ±10 per- 
cent due to all instrumental errors. The altitude resolution at 20 km was 
0.5 km for channel A and 0.7 km for channel B, values which are comparable to 
the 0.45 km value for the lidar profile. Details of the general operation and 
calibration of the particle counter have been described in references 5 and 13. 


TABLE I.- LIDAR AND DUSTSONDE MEASUREMENTS 
FOR TEXAS COMPARATIVE EXPERIMENT 


Date 

— 
Time, MDT 

Number of shots 

Lidar soundings 

May 7, 1974 

2035 

50 


2102 

50 


2223 

60 


2333 

300 

May 8, 1974 

0122 

70 


0158 

100 


0236 

130 


0328 

100 

Dustsonde launch 

May 8, 1974 

0200 



In figures 9 and 10 the aerosol backscatter ing-f unction profile is compared 
directly with particle counts N^ measured with the dustsonde. Error bars on 
f a represent 95-percent confidence intervals determined from the lidar signal 
statistics which vary according to the signal-to-noise ratio in each of the 
photomultiplier tubes. If the 1 -percent random error in the molecular density 
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Backscattering function 

Particle concentration (r > 0.15 (xm) 





were included, the total standard random error in f a would be 35 percent at 
20 km. Since initially only a relative profile comparison was sought, the 
dustsonde and lidar profiles in figures 9 and 10 were normalized at one alti- 
tude, 20 km. Inherent in this method of profile comparison is the assumption 
that the a a is constant with altitude in the stratosphere. If it is further 
assumed that the dustsonde and lidar are responsive to the same portion of the 
particle size distribution, then N a = N^. The comparison in figure 9 for 
particles with radii greater than 0.15 Mm shows that 1 particle per cm^ corre- 
sponds to a value of f a = 4 x 10“^ m _1 -sr _1 . Minima occur in both profiles 
between 15 and 17 km, and enhancements exist from 17 to 22 km. The strato- 
spheric aerosol layer exhibited a maximum concentration of 0.6 particle per cm^ 
at 19 km. If the standard error of 10 percent for is considered, the 

relation between the measurement techniques is not within the error bars above 
23 km. This discrepancy is addressed later in this section. Figure 10 shows 
the same comparison for particles with radii greater than 0.25 Mm. Here the 
correlation is better, especially at the higher altitudes, indicating that the 
lidar wavelength, 0.694 Mm, may exhibit sensitivity to the particle size dis- 
tribution with height. Particles detected in channel B exhibit correspondence 
when 1 particle per cm-* is equivalent to f a = 1 .8 x 10~8 m“^-sr“^. Values for 
O a derived from figures 9 and 10 are 4 x 10 -1 ^ m^-sr -1 and 1.8 x 10~ 14 m^-sr' 1 , 
respectively. At altitudes closer to the tropopause, these values of 0 a do 
not apply since cirrus particle composition, size, and shape are considerably 
different from those for the aerosols. 

Now, consider the two assumptions applied to the TCE data analysis: 

(1) that the 28 to 31 km region for normalization was aerosol free and (2) that 

the quantity 0 a was constant between 11 and 28 km. 

\ 

Effect of Normalization Criterion 

There is no evidence from figures 6 and 8 to suggest that a different alti- 
tude, between 10 and 31 km, would have been better than the 28 to 31 km range 
for normalization, although Russell et al. (ref. 3) present data that indicate 
an average underestimate in R of 3 percent if 30 km were chosen instead of 
12 km. Still, it is possible that no altitude between 10 and 31 km is really 
free of aerosol. 

Aerosol model calculations in reference 14 by Pinnick et al. (see their 
model N] with refractive index 1.40 - Oi) indicate that for N<j = 0.06 cm - -* at 
28.0 km in figure 9, one could expect f a = 5 x 10 - ^ m _ ^-sr - ^ (or R = 1.05) 
at that level. (Of course, this estimate assumes that the lidar and dustsonde 
are responsive to the same part of the aerosol size distribution; that is, 

Nd = N a .) By applying this same model at 1 1 km, where R is also a minimum 
and where the dustsonde data are more reliable, N<j = 0.78 cm - -* corresponds 
to f a = 4.5 x 10“9 m - l-sr -1 (or R = 1.04). A 5-percent increase in R 
at the normalization altitude results in a change in f a at 20 km of about 
3 x 10 -9 m _1 -sr -1 , an amount equal to 63 percent of the estimate of f a itself 
for a background aerosol layer. Since these aerosol model calculations yield a 
normalization uncertainty in R of about 5 percent at 30 km, this criterion 
may be the limiting systematic error in this present inference of absolute con- 
centrations of background stratospheric aerosols by lidar. 
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Figure 10.- Comparison of lidar and dustsonde profiles for particles 
with radii greater than 0.25 ym from TCE. 


Effect of Constant Backscatter Cross Section 


Data presented in figures 9 and 10 suggest that the normalized lidar mea- 
surements more closely approximate the larger particle profile (channel B) . 

The comparison with particles measured by both channels A and B is reasonable 
up to 23 km; above that altitude, the lidar return seems to underestimate the 
particle count in channel A. 


These results can be analyzed by evaluating the sensitivity of lidar back- 
scatter to various aerosol properties such as the refractive index and the 
effective cross-sectional area. In early 1974, Russell et al. (ref. 11) com- 
pared airborne impactor samples of prevolcanic aerosols to lidar returns at one 
altitude of 18 km. The refractive index n applied to the aerosol model in 
that comparison was 1.41 with no imaginary component, consistent with a concen- 
trated sulfuric acid composition. The particles were assumed to be spheres. 

The assumption of sphericity seems reasonable since Farlow et al. (ref. 15) 
report collections of only liquid aerosols between 12 and 32 km and for a range 



Figure 11.- Profile of aerosol number density ratio from dustsonde launch of 
May 8, 1974. (See eq. (A4).) 1 mb = 100 Pa. 
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of latitudes. By applying Mie theory to particles with refractive index 1.41, 
Russell et al. (ref. 11) calculated aerosol scattering functions which matched 
their lidar results at 18 km. For determination of aerosol profiles, the pos- 
sibility of a refractive index change must be considered in a a if the aerosol 
composition varies with height. Lazrus and Gandrud (ref. 16) reported aerosol 
measurements with height and latitude that showed the predominance of sulfate 
in the 1 2 to 28 km region. The bulk of this material is presumed to be concen- 
trated sulfuric acid. Therefore, the TCE scattering ratios above 11 km are 
assumed to have been affected most by variations due to changes in the particle 
size distribution rather than changes in composition. 

An altitude varying size distribution affects the effective a a , which 
relates backscattered lidar signals f a and particle number densities N,j. 

This relation is examined in appendix A, in which the slope of a power-law size 
distribution function is related to the aerosol number density ratios (channel A 
divided by channel B) measured by the dustsonde. Aerosol number ' density ratios 
for the TCE dustsonde are plotted in figure 11. At all altitudes above 10 km, 

the aerosol number density ratio is greater than 3.5; from 20 to 28 km, the 

ratio increases from 4 to 10. An increase of the aerosol number ratio with 
altitude above the main aerosol layer is normally observed (ref. 17) and can 
be explained by the theoretical arguments presented by Hamill et al. (ref. 12). 
The changing number ratio or size distribution causes the lidar backscatter ing 
efficiency or the effective f a (f a = 0 ^ 3 ) to decrease with altitude. By 
multiplying the appropriate efficiency factors y, derived in appendix A for 
the data in figure 11, with the curve of in figure 9, a profile of effec- 
tive f a is generated. That calculated profile is normalized at 20 km to the 

profile of measured f a of figure 9; both profiles are plotted in figure 12. 

The comparison from 23 to 28 km has been improved by as much as 30 percent. A 
similar calculation for data in figure 10 (r greater than 0.25 pm) yields an 
underestimate of f a above 23 km. Since no data exist on the exact size dis- 
tribution function for the particles measured during the TCE, the preceding 
calculations cannot be performed more accurately. Note that this final uncer- 
tainty due to cr a ( z ) affects both the scattering-ratio profile shape and mag- 
nitude, and it can vary with the laser wavelength (o a is a function of X). 

The effect of a variable aerosol scattering cross section can be important 
for volcanic as well as background stratospheric-aerosol layers (ref. 17). How- 
ever, McCormick et al. (ref. 4) have shown that the bulk of the stratospheric-™- 
aerosol backscatter is derived from the region below 23 km, where the aerosol 
number density ratios are normally less than 3.5. Thus, trends in integrated 
aerosol column densities for volcanic aerosol layers should not be affected by 
variations in a a with time. 


Background Aerosol Attenuation Profile 

Turbidity T is a term used to indicate the aerosol loading in the strato- 
sphere. It is defined as the ratio of the aerosol (B a ) and molecular (3m) total 
scatter functions, each with dimensions of inverse length; that is. 
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T 


(3) 


3a 

3m 



Figure 12.- Comparison between measured and calculated aerosol 
backscatter ing-f unction profiles. 


Elterman (refs. 9 and 10) infers turbidity values from his searchlight mea- 
surements of attenuation. It is also possible to derive a turbidity estimate 
and an attenuation profile from the TCE data. The relation, developed in 
app end ix B, between scattering ratio and turb idit y fo r a backgroun d a ero sol 
model is written as 


R 


1 .0 + 0.195T =1.0+— 

^m 


(4) 


As shown in appendix B, the coefficient of T in expression (4) depends on 
the aerosol composition and size distribution. A Deirmendjian Haze H model 
(ref. 18) of concentrated sulfuric acid was selected as representative of the 
stratospheric-aerosol layer in May 1974 because it was consistent with the 
background aerosol estimate of Harris and Rosen (ref. 19) and is equivalent to 
the power-law size distribution given in appendix A for those particles with 
radii greater than 0.2 ym. 
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Reported turbidity values of 1.4 at 20 km prior to the Mt. Agung eruption 
in the spring of 1963 and of 1.5 in November 1970 (ref. 20) correspond to scat- 
tering ratios of 1.27 and 1.29, respectively, for the aerosol model applied in 
expression (4). The derived scattering ratio for 1970 is somewhat higher than 
the lidar observations of Fox et al. (ref. 1) and others for 1970 and 1971 
(fig. 1), but this could be due partially to the presence of aerosols at 
normalization altitudes. If we assume that R = 1.05 (rather than 1.0) at 
28 to 31 km according to our discussion of the normalization uncertainty, then 
expression (4) applied to the adjusted TCE scattering ratio of 1.12 at 20 km 
yields a background turbidity of 0.62 for May 1974. Elterman (ref. 21) aver- 
aged 113 turbidity profiles obtained over New Mexico from September 1973 to 
July 1974, showing a value of T of 1.0 at 20 km. However, he reports an 
aerosol optical thickness between 12 and 25 km for April 1974 that is a factor 
of 1.3 less than his 10-month average. Therefore, a first-order estimate of 
his turbidity at 20 km at the time of the TCE is 0.78, or about 25 percent 
higher than the value inferred from lidar. This discrepancy is well within the 
absolute error estimates for both f a and the term 0.195T in expression (4). 

Scattering-ratio profiles can also be converted to yield aerosol attenua- 
tion profiles using the same aerosol model. Profiles of 3 a are obtained by 
combining data from expression (4) with a tabulation of the 3m profiles. If 
5 percent of the TCE lidar return signal at 28 to 31 km is assumed due to aero- 
sol (R = 1.05), as was estimated in an earlier section, then the curve f a in 
figure 9 is adjusted to give the aerosol .attenuation profile in figure 13. 

This profile is offered as an estimate of the aerosol extinction for back- 
ground stratospheric aerosols. 


CONCLUSIONS 

Just after the 1972 Wyoming comparative experiment (WCE) , a multiplexer 
was added to the Langley 48-inch lidar system which enabled the recording of 
a single stratospheric return (10 to 30 km) using three photomultiplier tubes 
to cover the signal dynamic range. Refinements in the calibration and data- 
analysis procedures were completed in order to minimize the occurrence of sys- 
tematic errors. The lidar system generated total stratospheric backscatter 
signals which when averaged yielded less than 2-percent random counting error_s 
for altitudes up to 30 km and over time periods of less than 1/2 hr. A lidar- 
dustsonde comparative experiment (TCE) was conducted in May 1974 at San Angelo, 
Texas, incorporating these updated system calibration and data-analysis pro- 
cedures. Uncertainties in scattering ratio R induced by the uncertainties in 
the molecular density profile and the attenuation model. were of the order of 0 
to 3 percent (average of 1 percent) and less than 1 percent, respectively. 

Both these uncertainties are assumed to be random. The root-mean-square 
average of all the random uncertainties is about 2.5 percent. For the TCE, 
the systematic uncertainty in R due to the clean-air normalization procedure 
is estimated to be about 5 percent. When these random and systematic uncer- 
tainties are applied to the TCE aerosol backscatter ing function f a at 20 km, 
they represent 35 and 63 percent of the value, respectively. 
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Comparisons between the dustsonde particle profile and the backscat- 

tering profile f a were within the measurement error bars for the two tech- 
niques below 23 km. The two-channel dustsonde data, however, indicate a chang- 
ing particle size distribution with height. This variation has been applied to 
the N <3 profile with the aid of a realistic aerosol model, and an effective 
f a profile was generated. As a result, the comparison between the two f a 
profiles is improved by up to 30 percent above 23 km. Since the bulk of the 
stratospheric-aerosol backscatter is derived from the region below 23 km, this 
variable size distribution effect is not very important in an analysis of trends 
in integrated aerosol column density. 

The TCE results indicate that individual measurements of scattering ratios 
above 1.025 can be significant provided that the random uncertainties due to 
noncoincident molecular-density profiles are avoided. Absolute values of the 
aerosol backscatter function f a are further affected by uncertainties in the 
profile normalization process, and this error must be considered when relative 
changes (seasonally or yearly variations) in aerosol backscatter are evaluated. 
Attempts to extend the altitudes from which lidar data are retrieved to greater 
than 30 km in order to decrease the normalization error are not advised unless 
methods for obtaining accurate molecular-density profiles from those same levels 
are also available. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 11, 1979 


APPENDIX A 


AEROSOL CROSS-SECTION CALCULATIONS 

In order to test the sensitivity of lidar returns to the size distribu- 
tion, assume a power-law size distribution 


ON 

— = Cr - < v+1 ) (Al) 

dr 

where N- is the number of particles per cm-*, r is particle radius in micro- 
meters, and C is a coefficient proportional to the total number of particles 
in the size distribution. (A power law is reasonable for 0.1 to 1.0 Mm radius 
particles.) The slope of the distribution is given by the exponent v. Val- 
ues (in m?-sr -1 ) of the backscatter ing cross section o a are given by the 
expression 


dn ^a 
dft N a 


(A2). 


where ft is a solid angle and u is the geometric cross section. Values for 
a a are derived from Mie theory as a function of v and of size parameter a 
where 


a 


2irr 

T~ 


(A3) 


Values for a a have been tabulated by McCormick et al. (ref. 22) for laser 
wavelengths.- Numerical integrations to obtain dn/dft over the size- 
distribution equation (Al) are performed for a wavelength of 0.694 ym, a 
refractive index of 1.41, and particle radii limits of r-j = 0.15 and 0.25 Mm 
to t 2 - 3.0 ym. Sample tabulations are given in table Al. 


The size distribution slope v can be related to the ratio of the number 
of particles in channel A, represented by Ng.iS/ to those in channel B, repre- 
sented by Nq 25» as follows: 


N 0 . 1 5 
N 0.25 


✓*00 

*A).i 


r" ( v+1 ) dr 


5 ym 


I 


r-( v+ Udr 


0.25 ym 


(A4) 


The left-hand side of expression (A4) is determined from measurements of N^ by 
the University of Wyoming dustsonde, and selected results appear in table All. 
Expression (A4) is considered valid, due to the dustsonde measurement technique. 
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N 0. 15/^0.25 


m 


?ERENTIAL CROSS SECTION 
TOR BACKSCATTER 


dri/dft, 

m^-sr 

-1 

r 2 = 3 - 

0 ym 


1 .52 x 

10 _1 

4 . 

7.77 x 

10 -1 

5 

4.50 



2.93 



r 2 = 3. 

0 ym 


5.18 x 

10 _1 

4 

3 . 25 



2.23 



1 .66 




:kscatter on size distribution 
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APPENDIX A 


for particle ratios greater than 3.0 (ref. 13). A lidar aerosol return is 
determined by the product of the relative backscatter function per particle 
dri/dft and the number of particles N a . From table AI, one can then form the 
ratio of that product for each of the two size ranges measured with the dust- 
sonde. As an example, v = 3.0 corresponds to N 0 . 15 /N 0.25 = 4 - 65 and 



dn\ 

d W0.25 


0.24 


(A5) 


The product of the two ratios, equations (A4) and (A5) , yields the magnitude of 
the lidar signal return to be expected from particles in channel A relative to 
that expected from channel B. The inverse of that product is the relative 
efficiency for backscatter, y = fo.25/ f 0.15» for fche s * ze distributions ratios. 
When y is less than 0.5, the ruby lidar receives more signal from the popula- 
tion of particles between 0.15 and 0.25 ym than from the population larger than 
0.25 ym. 

The power-law size distribution function does not apply so well to parti- 
cles less than 0.15 ym in radius. Some particle measurement techniques even 
indicate decreasing or constant particle numbers below about 0.05 ym in radius. 
Calculations for more realistic distributions of the particles smaller than 
0.15 ym indicate that only a very small error is incurred in the above analysis 
by neglecting them altogether . 
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APPENDIX B 


TURBIDITY CALCULATIONS 

A measure of atmospheric particulate loading is the turbidity T, defined 
as the ratio of the aerosol and molecular total scatter functions 


e a 



Turbidity measurements are presented for a wavelength of 0.55 ym in refer- 
ence 20. They are converted to 0.694 ym to correspond with the lidar results 
of the present study by 

/ \ 2 * 5 

I 0.694 \ 

T 0 . 694 = T 0 . 55I - -- J = 1 -77T 0>5 5 (B2) 

where 8 m varies according to X -4 and $ a varies according to X -1 -5. 

The aerosol wavelength dependence is taken from figure 7 of Pinnick et al. 

(ref. 14). A backscatter turbidity T(180) is measured by lidar and is 
defined as 

P ( 1 80) 

T(180) 0.694 = " T 0. 694 • ( B3 ) 

I • b 

where 1.5 is the Rayleigh scatter phase function. The aerosol backscatter ing 
phase function P(180) is given by Elterman et al. (ref. 20) to be 0.232 for 
a silicate Deirmendjian Haze H aerosol of refractive index 1.54 (ref. 18), 
while Russell et al. (ref. 11) use P ( 1 80) = 0.165 for a sulfuric acid Haze H 
model (refractive index of 1.41). Since sulfuric acid is the preferred compo- 
sition, application of this latter value to equations (B2) and (B3) yields 

T(180) o, 694 = 0.195T 0 .55 (B4) 

The expression for the scattering ratio can be written as 

R = 1.0 + T( 180) 0 . 6 94 = 1.0 + 0.195T 0-5 5 (B5) 

In table BI the sensitivity of the relation between Elterman' s turbidity 
and the scattering ratio is presented by considering three different haze mod- 
els. Since the coefficient of turbidity varies sharply with the chosen strato- 
spheric model, expression (B5) should only be applied to a relatively unper- 
turbed, nonvolcanic aerosol layer. 
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